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Over the last years, superconductivity has been discovered in several families of iron-based com- 



;_( pounds. Despite intense research, even basic electronic properties of these materials, such as Fermi 

-(— > 

. surfaces, effective electron masses, or orbital characters are still subject to debate. Here, we address 

-(— > 
J5 an issue that has not been considered before, namely the consequences of dynamical screening of the 

H 

I Coulomb interactions among Fe-d electrons. We demonstrate its importance not only for correlation 

f— satellites seen in photoemission spectroscopy, but also for the low-energy electronic structure. From 

r^ our analysis of the normal phase of BaFe2As2 emerges the picture of a strongly correlated compound 

with strongly doping- and temperature-dependent properties. In the hole overdoped regime, an inco- 

K^ herent metal is found, while Fermi-liquid behavior is recovered in the undoped compound. At optimal 

"^^ doping, the self-energy exhibits an unusual square-root energy dependence which leads to strong band 

^-H renormalizations near the Fermi level. 

I • Most known superconductors can be attributed to one of two classes: The first class comprises materials such 

^D as MgB2, classified as "weakly correlated" in the sense that one-electron theories work well for the description of 
T-H the basic electronic properties, and the superconducting pairing mechanism has eventually been understood in this 

^ framework. Materials of the second class exhibit sometimes spectacular failures of the one-electron picture, leaving 
/\ unclear even the theoretical language in which a theory of the pairing mechanism should be formulated. High-Tc 
^ cuprate superconductors fall in the latter category. The parent compounds of these materials are Mott insulators 
due to strong Coulomb interactions that localize the electrons on the Cu sites. The role of electronic correlations in 
the recently discovered iron-based high-Tc superconductors 1 is less clear, and apparently dependent on the specific 
family, as well as on doping, substitutions or pressure. In this work, we address the prototypical compound of the 
so-called "122-family" , BaFe2As2, which exhibits superconductivity under pressure |21I3] or hole- as well as electron- 
doping [U \S\ . Various experimental probes - angle-resolved and angle-integrated photoemission spectroscopy [SHU] , 
optics and transport, Raman spectroscopy, neutron experiments, NMR - have been employed to characterise the 



electronic properties [T3]- The Fermi surface consists of two concentric hole pockets around the F point, and elliptic 
electron pockets around M 14J. Experimental estimates of the (doping-dependent) mass enhancements vary from 
about 1.4 [15 to 5, at least for the orbital pointing towards the As-sites [IS]. The orbital character of these pockets 
are still subject to debate, but there seems to emerge a consensus about stronger correlation effects for holes than 
electrons. 

Understanding the low energy electronic structure, in particular Fermi surface nesting, orbital character and mass 
enhancements, is a prerequisite for assessing possible pairing mechanisms JT51 [TTHT^ . While the field of electronic 
structure calculations for correlated materials has made tremendous progress in recent years, most notably thanks 
to the combination of electronic structure and many-body techniques, the discovery of iron pnictides provides a 
serious challenge. The combined "LDA+DMFT" scheme builds on density functional theory within the local density 
approximation (LDA) to construct realistic many-body Hamiltonians, which are solved using dynamical mean field 
theory (DMFT) [5Dl - f22] . While a number of interesting applications of this scheme to iron pnictides (see e.g. [71 [2^[M| ) 
have emerged, it has become clear that one of the bottlenecks is the determination of the Coulomb matrix elements 
("Hubbard parameters" U and J), that parametrize the energetic cost associated with the distribution of electrons 
among the localized Fe-d orbitals. Recent works try to extract these values from constrained LDA calculations 
[23] . GW-inspired methods [IB", or from random-phase approximation (RPA) based schemes. The constrained RPA 
approach T] takes into account the dynamical nature of screening, so that the "Hubbard f7" for the Fe-d states becomes 
a frequency-dependent object. None of the theoretical studies so far have treated dynamical Coulomb interactions 
at the level of the actual many-body calculation. Here, we demonstrate the importance of the frequency-dependence 
for the low-energy electronic structure, mass enhancements and lifetimes, as well as for a description of correlation 
satellites seen in photoemission spectroscopy. 
Partially screened Coulomb interaction 

The constrained RPA result for the average intra-orbital Coulomb repulsion U of BaFe2As2 is shown in Fig. [T^. 
The real part ranges from the static value C/q = ReC/(a; = 0) = 3.6 eV to the bare interaction V of about 20 eV 
at large w, and the frequency dependence resembles typical U{io) in transition metals []J|6j, where there is no well- 
defined plasmon excitation. Instead, the main high-energy excitation with a peak in lmU{uj) at ^ 26 eV forms 
a broad structure extending down to a few eV, implying that the plasmon excitation overlaps strongly with the 
one-particle excitations, a reflection of the semi-itinerant nature of the Fe-Sd electrons. U{io) represents a partially 
screened Coulomb interaction for the Fe-d states, which accounts for screening by all degrees of freedom except the 
Fe-d states themselves. Fourier transformation yields a Hubbard interaction term j ^ drdr' N{t)U{t ~ t')N{t'), 
with U{t) = V5{t) -\- t/rctardod(''"), which dcscribcs the local Coulomb interaction including retardation effects due to 
quantum fluctuations that screen the bare V (see Fig. fit and Methods Section). 

In a standard LDA+DMFT calculation without dynamical effects, the interaction Uq would result in a rather 
weakly correlated picture. This is demonstrated in Fig. [2J3, where the Fe-d spectral function from such a static-[/ 
LDA-I-DMFT calculation is shown. Interaction effects lead to a moderate renormalization of the Fe-d states, with a 
mass enhancement of 1.6. A comparison to the LDA density of states in panel (a) shows that the peaks at —3 eV and 
1 eV are weakly renormalized band states. No Hubbard satellites or other correlation features appear, in agreement 
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FIG. 1: Frequency dependent interaction for BaFe2As2 from constrained RPA. Panel (a) shows the real and imaginary parts of 
the average intra-orbital Coulomb interaction as a function of frequency. V is the unscreened (bare) Coulomb interaction and 
Uq the static (partially screened) interaction. \vs\U (iS) j uP' , plotted as a red line, features prominent peaks at 3.8 eV, 6.1 eV, 
16 eV, and 26 eV, and smaller humps at 10 eV and 12 eV. Panel (b) shows the spectral function of the bosonic propagator B. 
The peaks of this function, inherited from the structures in Im(7(cj)/ctj^, determine the energies of satellite structures in the 
d-electron density of states. A schematic plot illustrating the difference between a static-L'^ (yellow) and a dynamic-?7 (blue) 
calculation is shown in panel (c). In a static-[/ calculation, up-spin and down-spin electrons interact locally in time with the 
partially screened Coulomb interaction f/o- In the dynamic-f/ calculation, the instantaneous interaction is given by the bare 
V , while screening effects lead to an attractive retarded interaction C/i.otardcd(''")- 



with previous studies [231 HI] • 

The new aspect of our work, compared to previous simulations, is the treatment of the full frequency dependence of 
the interaction. In Refs. [T5l ll4l it was shown in the context of simple model calculations, how a frequency dependent 
interaction can be incorporated efficiently into Quantum Monte Carlo simulations within the hybridization expansion 
impurity solver scheme J15j . Here, we generalize this technique to multi-orbital systems and realistic materials. The 
effect of t/(w) is to dress the Fermionic propagators with a bosonic propagator 



B(t) = exp \-K{t)\ , 

where K(t) is the twice integrated retarded interaction. In terms of ImC/(w) and a factor 6(r, w) — cosh [(r 
^)aj]/sinh \pf\ with bosonic symmetry, we can write |14j 

[°° dtu ImC/(w) 



(1) 



K{r) 



[6(r,^)- 5(0,6.)]. 
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It is evident from Eq. (11 1 that, as far as structures in the frequency dependent interaction are concerned, the relevant 
function to analyze is '°^2 . This function is plotted as red dash-dotted line in Fig. ik. Besides a first peak at 3.8 eV, 
which comes from the rapid decay of Im[/(w) at small frequencies, there are prominent peaks at 6.1 eV, 16 eV and 
26 eV, as well as smaller features at 10 eV and 12 eV. We have traced the origin of the 6.1 eV feature to transitions 
from occupied d-states to states in the energy window [6 eV:7 eV], which have predominantly Ba character. The 12 
eV feature is due to transitions from As-s semicore states to empty Fe-d states close to the Fermi level. 
Density of states and high-energy satelHtes 




FIG. 2: Effect of (7(cj) on tlie d-electron spectral function. Panel (a) shows the total LDA density of states for BaFe2As2 
with the p- and d-electron contribution shaded in red and blue, respectively. Panel (b) shows the d-electron spectral function 
obtained from the DMFT calculation at T = 290 K and optimal hole doping (x — 0.2 per Fe), using the static value Uo (red 
curve) and the full frequency dependent U{ijj) (black curve). The inset shows the high-energy tail of the d-electron spectral 
function. Arrows indicate the position of satellites predicted by the dynamic-t/ calculation. In panel (c), we plot the total 
spectral function (p- and d-contribution) of optimally doped BaFe2As2 obtained from DMFT. 



An interesting question is how the structures in U{uj) affect the properties of the material. As discussed in Ref. \T^ 
and in the Supplementary Material, the Green function for a system with frequency dependent interactions U{u}) 
factorizes in the atomic limit into the atomic Green function for the system with static interactions Uq and a bosonic 
factor, which is precisely the propagator B defined in Eqs. (fl]) and 111]). Inspired by this atomic limit property, we 
introduce the auxiliary Green function GauxiUary by G{t) — G aMyS&Bxyij^^ij^ and compute the spectral function of the 
system as the convolution of the spectral functions of the auxiliary propagator, Paux(w), and of the spectral function 
of the bosonic propagator, ps(a;) (see Eq. (9) in the Supplementary Material). 

For systems with a single well-defined plasma frequency the effect of this convolution would be to replicate the 
low-energy spectral features of Paux(i^) (with exponentially decreasing weights), displaced by multiples of the plasma 
frequency. These satellites correspond to processes where in addition to the one-electron addition or removal process a 
certain number of plasmons are emitted or absorbed. In the present case, the bosonic spectrum, plotted in Fig. [Tb, is 
more complex than just a single plasmon delta-function, since it inherits the structures of linU (uj) / lo'^ . Nevertheless, 
sharp features present in pb{^) lead to replications of the structure of the low-energy spectral density. 

The low-energy part of the d-electron spectral function from the dynamic-C/ calculation is shown in Fig.[2]D, and the 
total spectral function {p- and d-contribution) in Fig. [2]:. In the dynamic-C/ calculation, spectral weight is shifted to 
high energies, which leads to a reduction in weight at low energies, compared to the static-C/ result. More importantly, 
we see that the structures arising from p-d hybridization are shifted closer to the Fermi energy, and the peak near 
a; == is strongly renormalized. Thus, the explicit treatment of the strong Coulomb repulsion at large frequencies 
has a substantial effect, even on the low-energy properties of the system. The large quantitative effect of U{io) is 
remarkable given the fact that the real part of the frequency-dependent interaction (Fig. fTk) remains of the order of 
Uq or lower up to w ~ 15 eV. 



The sharp low-energy peak in the d-electron spectral function results in weak, but well-defined satellites, as discussed 
above. The inset of Fig. [2b shows the high energy tail of the occupied part of the spectrum, with arrows marking 
the most prominent satellites at —6.1 eV, —12 eV and —16 eV. This physics is of course not contained in a static-C/ 
calculation or any other previous theoretical work on pnictides. The observation of satellites at —6.5 eV and —12 eV 
was emphasized in the photoemission study of Ding and collaborators [S]. While the feature at —12 cV was argued 
in Ref . 1331 to be an As-4s line, the latter work confirms a hump in the d-electron spectral function around —6.5 eV. 
Our calculation suggests that a d-feature, originating from the structure in the frequency dependent interaction, is 
superimposed to the As-4s spectral contribution. The —16 eV feature is probably not visible in experiments, because 
it overlaps with Ba-5p states, while a satellite which we predict at —3.8 eV is masked by structures arising from p-d 
hybridization. 
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FIG. 3: Non-Fermi liquid behavior of the self-energy. Panel (a) shows the imaginary part of the self-energy on the Matsubara 
axis for optimal hole doping (x — 0.2 per Fe), both for the dynamic-[/ calculation (red circles) and the static-17 calculation (blue 
diamonds). Solid and dashed lines are fits of the low-frequency behavior to the function — ImE(ia;„) = A{iJn)°' ■ The inset shows 
the low-frequency behavior of the dynamic-f/ result for different temperatures. As the temperature is raised, the extrapolation 
to iOn ^ yields a non-zero intercept, which indicates that even excitations at the Fermi level exhibit a finite lifetime. Panel 
(b) shows the low energy behavior of the self-energy as a function of doping. Fermi-liquid behavior is found in the undoped 
and electron-doped compounds, while a non-zero intercept appears in the overdoped case. Panel (c) shows a sketch of the 
phase diagram in the space of temperature and doping. The blue region indicates Fermi-liquid behavior, while yellow indicates 
a frequency-dependence of the self-energy which is not compatible with Fermi-liquid theory. The left dashed line marks the 
boundary of the crossover region, where the exponent a starts to deviate from 1. The right dotted line corresponds to q = 0.5, 
which marks the "spin freezing" transition. To the right of this line, an incoherent metal phase with non-zero intercept of 
ImE is found. The experimentally measured phase-diagram with superconducting (SC) and spin-density-wave (SDW) ordered 
phases is indicated by black dotted lines. Full substitution (KFe2As2) corresponds to a; = 0.5. 



Non- Fermi liquid properties of the metallic phase 

In a Fermi liquid, the imaginary part of the Matsubara axis self-energy exhibits a linear regime at low energy, whose 
slope is directly related to the quasi-particle mass enhancement. However, as shown in Fig. [3^, for optimally doped 
BaFe2As2 (x = 0.2 hole doping per Fe), and at the temperatures of our simulations, we do not observe this behavior: 
the self-energy behaves as — Iml](i(u;„) — A{iuJn)°' with a w 0.56 at T = 145 K, while a frequency dependence of the 



form — ImS(iaj„) ^ B + A{iujn)°', with non-zero intercept (scattering rate) B is found at T = 290 K and T = 580 K. 

A similar phenomenon has recently been observed in the metallic phase of a 3-orbital model 35J, and has been 
dubbed "spin-freezing" transition, because the intercept B is due to scattering off static (but disordered) local mo- 
ments. The exponent value a — 0.5 (and simultaneous onset of an intercept B > 0) was found to mark the transition 
into the spin frozen, or incoherent metal regime. Large deviations from Fermi-liquid behavior thus appear in the 
vicinity of this transition. The phenomena described in Ref. 1351 seem to be a generic property of the metallic phase in 
multi-orbital systems with large Hund's coupling. Our results in Fig. [3^ indicate that optimally doped BaFe2As2, at 
temperatures somewhat above Tc, is close to the spin-freezing transition, and that increasing the temperature shifts 
the material across the spin-freezing boundary into the incoherent metal regime. 

A square-root like rather than linear low-frequency behavior of ImE means that Landau quasi-particles and effective 
masses cannot be properly defined. The square root implies that the bands very close the the Fermi level are much 
more strongly renormalized than those further away, consistent with the behavior seen in Fig. [2] This non-Fermi 
liquid (NFL) property may explain why different photoemission experiments arrive at considerably different estimates 
for band rcnormalizations, and why it has been so difficult to reach a consensus on the importance of correlations in 
BaFe2As2 (and, possibly, other pnictides). 

While Fermi-liquid properties may eventually be recovered in simulations at low enough temperature (see Fig. ^) , 
such behavior is cut off in real materials by the onset of the spin-density-wave or superconducting phase. Thus the 
NFL and incoherent metal regime dominates the physics in the whole temperature range of relevance to this study. 
We also note that NFL behavior has been proposed for LaFeAsO and FeSe, based on static-C/ calculations PUIITTIIST] . 

Our result from the static-C/ calculation is plotted in Fig. |3^. We see that taking into account the frequency 
dependence leads to a substantially enhanced self-energy, and a much larger frequency range over which the NFL 
exponent is valid. Fitting the self-energy from the static calculation to Iml](zw„) = A{iiOn)°' yields a ~ 0.75, which is 
much closer to Fermi-liquid behavior. In this sense, the dynamic-[/ calculation, by increasing the interaction effects, 
shifts the material closer to the spin freezing transition, and the proximity to this transition line results in a sensitive 
dependence of simulation results on parameters such as temperature and doping level. 

Figure [3]d shows the low-energy part of the self-energy of BaFe2As2 at T = 145 K for different dopings. Increasing 
the hole doping leads to a large scattering rate while reducing the doping reduces the slope of — ImS and hence the 
band renormalization. In fact, for undoped and electron doped materials we find a Fermi-liquid behavior at the lowest 
temperatures, with mass renormalization factors 2.6 (undoped) and 1.8 (electron doping x — 0.2 per Fe), respectively. 

Our results on the NFL behavior are summarized in Fig. [3J;, which sketches the phase-diagram in the space of 
temperature and doping. The incoherent metal regime is shown in yellow, and the Fermi-liquid region in blue. The 
dotted line is defined by the power-law exponent a — 0.5 (with simultaneous onset of static magnetic moments) and 
represents the boundary of the spin-frozen region in the sense of Ref. [351 The dashed line marks the temperature and 
doping level below which Fermi-liquid behavior is recovered. 
Momentum-resolved spectral functions and comparison to experiment 

Spectral functions for optimally doped BaFe2As2 are shown in Fig. |4j The orbitally resolved d-electron spectral 
functions along the path F-X-M-F are plotted in a narrow energy range in panel (a). A comparison with the original 
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FIG. 4: Orbitally resolved spectral functions for optimally doped BaFe2 As2 at T = 145 K. Panel (a) shows the spectral functions 
along the path F-X-M-F and compares the DMFT result to LDA calculations (dashed). The bands are strongly renormalized 
near the Fermi energy. Panel (b) shows the total p- and d-contribution to the spectral function in a wide energy window. 
The p-states hybridize with the d-states mainly in the unoccupied part of the spectrum. This is because it is precisely this 
hybridization that pushes the d-states up in energy. Panel (c) shows a cut through the hole pocket near F, for different values 
of fez, and resolved into x^ — y^ and xy + yz contributions. Dotted lines indicate the LDA bands. 



LDA band structure (dotted lines) shows the strong band renormahzation at low energies, a consequence of the 
"square root" behavior of the self-energy near the spin freezing transition. We also see that the z^ orbital has no 
weight at the Fermi energy, which at first sight might suggest a simpler model for BaFe2As2, involving less d-orbitals. 
However, due to the sensitive doping dependence of the spin-freezing phenomenon, such models may have properties 
which are very different from those of our five-orbital model. The spin freezing transition produces dramatic effects 
in the metallic phase away from half-filling^ where it leads to strong band renormalizations, even though the material 
is not close to any Mott insulating phase. In the five-orbital description, optimally doped BaFe2As2, with ^ 5.8 
electrons per Fe, falls into this NFL region (far away from the 6-electron Mott phase). On the other hand, in a nearly 
half-filled four-orbital model, NFL effects comparable to those demonstrated in Fig. [3] are expected to appear only in 
the vicinity of the Mott insulating phase. 

Panel (b) plots the total p- and d-contribution to the spectral function along the same path. Remarkably, the p- 



bands (energy range [—5 eV:— 2 eV]) overlap with the d-bands (energy range [—2 eV:2 eV]) mainly in the unoccupied 
part of the spectrum. This property, which can already be seen at the LDA level [S^, is amplified in the DMFT 
description, since the empty d-states are much better defined than the filled states whose spectra are smeared out by 
the self-energy effects described above. It may lead to an additional asymmetry between electron and hole doping, 
beyond the strong doping asymmetry implied by the spin-freezing phenomenon. 

A close-up view of the hole-pocket near F is shown for several values of kz in Fig. |4];. Dashed lines again indicate 
the LDA bands. While the renormalization of these bands is very large, the fc^-dispersion of the Fermi surface is 
found to be weak in the DMFT description. The outer Fermi surface has predominantly x'^ — y^ character and the 
inner Fermi surface xz + yz character. 

A comparison of low-energy momentum-resolved spectra to photoemission data is presented in Fig. pj Panels (a) 
and (b) show the spectral function along the path F-X-M-F, for optimal hole doping x — 0.2, T — 145 K (top) and 
T = 290 K (bottom). Angle resolved photoemission data from Ref. [5]are indicated by green dots. (The photoemission 
spectra were measured at T = 50 K and T — 145 K, respectively). The agreement near the F, X and M points is 
remarkably good, which shows that the very strong low-energy renormalization implied by the NFL self-energy is 
consistent with experiment. Increasing the temperature leads to a smearing of the bands, a result of the temperature 
dependent intercept of ImS shown in the inset of Fig. [sk. Panels (c) and (d) illustrate the strong doping dependence 
of the spectra. In panel (c), we show the result for the undoped compound at T = 145 K, which is in the Fermi- liquid 
region of the schematic phase diagram (Fig. Isb) . Indeed, our calculation produces well-defined bands with a modest 
renormalization consistent with the factor 2.6 extracted from the slope of the Matsubara axis self-energy. In panel 
(d), we plot the spectral function of the overdoped sample {x — 0.4), which falls into the incoherent metal phase 
characterized by a non-zero intercept of ImS. As a result, the bands are smeared out. We also note that the electron 
pocket at M has disappeared. These theoretical predictions are consistent with measurements on the end member of 
the series, KFe2As2, reported in Refs. HOI BTl Panels (e) and (f) of Fig. [5] compare the low-energy band structure in 
the undoped compound at T = 145 K to photoemission data taken from Ref. [HI (green dots) and the LDA bands 
(black lines). In this case, our calculation underestimates the band renormalization somewhat. However, given the 
large differences in the band renormalizations between the optimally doped and undoped compounds, we consider 
the agreement with the theoretical calculation satisfactory. We note that the discrepancy with the experiment is less 
than 0.05 cV, and thus much smaller than in any previous theoretical work. 

The spin-spin correlation function is known to exhibit an unusually slow (imaginary time) decay near the spin 
freezing-transition, while there is no particular anomaly in the orbital correlation functions [55] . Our observation that 
the maximum Tc in BaFe2As2 is reached in the vicinity of the spin-freezing transition may therefore suggest that the 
spin degrees of freedom, rather than the orbital degrees of freedom play an important role in the mechanism leading 
to high temperature superconductivity. 
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FIG. 5; Doping and temperature dependence of the low-energy spectral functions, and comparison to photoemission exper- 
iments. Panels (c) and (f) show the result for optimal doping {x = 0.2) and compare the renormalized bandstructure to 
photoemission data taken from Ref. 13 We plot the result for two temperatures, T = 145 K (top) and T — 290 K (bottom), 
to illustrate the smearing of the bands at elevated temperature. Panels (b) and (e) illustrate the doping dependence. The 
undoped compound is in the moderately correlated Fermi-liquid regime, while the overdoped compound (x = 0.4) is in the 
incoherent metal regime, where the scattering off static moments washes out the bands. Panels (a) and (d) compare simulation 
results for the undoped compound at T = 145 K to photoemission data from Ref. [IGJ and the LDA band structure (lines) . 

METHODS 



Our scheme can be viewed as an extension of the combined LDA+DMFT method |21l I22j to dynamical interactions, or 
as promoting extended dynamical mean field theory (EDMFT) j42l I43j to a realistic scheme by combining it with input from 
electronic structure calculations within the LDA and constrained RPA frameworks. Alternatively, our method can be considered 
as an approximation to a full GW-I-DMFT calculation [44|. The present scheme simplifies the general formulation in so far 
as only a local self-energy is calculated (as in DMFT) and two-particle quantities are calculated at a non-selfconsistent level. 
Thereby practical calculations become feasible even for complex multiband systems such as the iron pnictide compounds. 



10 
In a Hamiltonian formulation, we can write the multi-orbital model with dynamically screened interactions as 
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+ ^ dij [Ai„(&]"„ + bi^)Y^ riima + i^hl^h^'^ . (3) 

i ma 

It consists of a multiorbital Hubbard-Hamiltonian of the form usually treated within LDA+DMFT '22] , albeit with the bare 
Coulomb interaction V entering the Hubbard-terms, and an additional bosonic Hamiltonian with bosonic modes coupling to 
the total electronic occupations of the atomic sites. Here, a\^^ creates an electron on atom i in orbital m with spin a and 
nimtj is the number operator for such electrons. For computational reasons, we have restriced ourselves to density-density 
interactions. The bosonic part of the Hamiltonian describes the coupling of the electronic degrees of freedom (via the total 
charge Ni = X^mcr "'»™o" '-'^ ^^^^ * ^'^'^ some coupling constant Xiio) to bosonic modes. These bosonic modes represent both, 
collective plasmon oscillations and modes parametrizing one-particle screening processes. 

We solve the multi-orbital lattice problem using dynamical mean field theory (DMFT) [20], which maps it to a self-consistent 
solution of a five-orbital quantum impurity model. This local approximation, and the integration over the bosonic degrees of 
freedom, leads to an action of the form 

a{r) 



+ 2J 


/ dr 


E Vmm'nma{T)r 


1 

+ 2J 


( dr 


^corrcl. orb.) 

/ ^ I, ^mm' ~" 'J mm 
(corrcl. orb.} 



(corrcl. orb.) 

+ ^ I drdr'7V(r)f/,etardcd(T - T')iV(r'), (4) 

with — Imf7i.otardcd(^) = 7i"Aj , aud the partially screened interaction f/o = V^ — 2 J doj-^. 

The frequency dependent (or retarded) U can be viewed as a systematically downfolded interaction, stemming from a static 
Hamiltonian including all - even high-energy - degrees of freedom, in the sense of Ref. [T| In practice, we calculate it from 
the constrained RPA method as described in the Supplementary Material. The impurity model (with frequency dependent 
interactions) is solved using a Monte Carlo method based on a stochastic expansion of the partition function in the impurity-bath 
hybridization [T4l[l5]. Details can be found in the Supplementary Material. 
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Supplementary Information 

These supplementary notes describe the methods used in our ab-initio simulation of BaFe2As2 with dynamically 
screened interactions. 

Frequency dependent interaction from constrained RPA 

In the constrained RPA (cRPA) method [3 [5] we first define a one-particle subspace {ipd} of the low-energy space, 
which we call the "d subspace", and label the rest of the Hilbert space by {V'r} ("?" subspace"). In the present case, 
we choose 10 states having strong Fe-3d character as the d subspace (see Fig. [6]). We define Pd{r , r' ; uj) to be the 




FIG. 6: Electronic band structure of BaFe2As2. The black lines show the LDA band structure, whereas the red dots are 
interpolated bands in the d subspace obtained using the maximally localized Wannier function procedure. 



polarization within the d subspace and P(r, r'; uj) as the total polarization. The rest of the polarization P,. = P — Pj, 
is not the same as the polarization of the r subspace alone because it contains polarization arising from transitions 
between the d and r subspaces. The physical idea behind the cRPA method is that the Hubbard U defined as an 
effective interaction of the electrons in the low-energy effective theory must be such that when it is screened by the 
polarization of the low-energy states, it should be equal to the fully screened interaction W of the whole system. We 
thus define the partially screened Coulomb interaction Wr{r,r';uj) by 



Wr{uj) = [l-VPr{uj)]-^V, 



(5) 



with V denoting the bare Coulomb interaction. The d states are entangled with others in the present system. In 
order to distinguish between Pd and Pr, we apply a disentangling procedure f3^. 

The Hamiltonian for our DMFT calculations contains ten bands around the Fermi level with dominantly Fe-d 
character and also six bands (located at —6 to —2 eV) mainly coming from the As-p bands. The Fe-d orbitals are 
treated as correlated, whereas the As-p orbitals are assumed to be non-interacting. The Hubbard U interactions for 
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the Fe-d orbitals are defined as matrix elements of Wr in the maximally localized Wannier function (MLWF) [U [S] 
basis for the Fe-d orbitals. More technical details on the estimation of U can be found in Refs.[6]and[7l 




20 30 

co[eV] 

FIG. 7: Frequency dependence of the fully screened Coulomb interaction W{io) for BaFe2As2. The dashed lines show the 
result for the partially screened Coulomb interaction U{uj). 



The frequency dependence of the partially screened interaction U{uj) is plotted in Fig. la of the main text. To 
illustrate the difference to the fully screened interaction W{uj) = [1 — VP{uj)]^^V, we plot W{uj) in Fig. ^ Note that 
R.eU{uj = 0) = 3.61 eV, while ReT4^(aj = 0) = 0.975 eV. The Hund coupling parameter J is approximately 0.8 eV and 
shows little frequency dependence. 

The frequency dependence of the Hubbard U reflects retardation effects arising from the r subspace that has been 
eliminated in the low-energy effective action. The derived action contains p and d orbitals with the onsite U{uj) term 



on the d orbitals. In the d-manifold with orbitals (xy, yz, z^, zx, x^ 
same and opposite spin are 



?7r„(J = 0.8) 



{/r„(J = 0.8) 



/3.6I 2.38 2.20 2.38 2.93 

2.38 3.61 2.74 2.38 2.38 

2.20 2.74 3.61 2.74 2.20 

2.38 2.38 2.74 3.61 2.38 

2.93 2.38 2.20 2.38 3.61 

0.00 1.76 1.49 1.76 2.58 ^ 

1.76 0.00 2.31 1.76 1.76 

1.49 2.31 0.00 2.31 1.49 

1.76 1.76 2.31 0.00 1.76 

2.58 1.76 1.49 1.76 0.00 



j/^), the static values of the interaction for 



(6) 



(7) 



and we take the frequency dependence shown in Fig. la of the main text for all components. In Eqs. ^m and m\ we 
have symmetrized the matrix elements obtained from cRPA, by computing the Slater paramters i*"", F^ and F"^ from 
the averaged C/"""" — 2.70 eV and J = 0.8 eV, and reconstructing the intra- and inter-orbital interactions using the 
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procedure described in Ref. |H1 This symmetrization procedure allows us to use an orbital independent double-counting 
term. 

Approximating the Coulomb interaction by the density-density terms leads to an underestimation of the spin 
fluctuations or, equivalently, to an overestimation of the role of the Hund coupling J. We found, in agreement with 
previous studies on other iron pnictides |10L lllj that the simulation results depend very sensitively on the choice of 
the parameter J . To obtain meaningful results for the low-energy spectral function, we thus reduced the value of the 
Hund coupling to J = 0.675, which (using tj'^'^ ~ 2.70 eV) gives the following interaction matrices: 

/^3.61 2.57 2.42 2.57 3.03^ 
2.57 3.61 2.88 2.57 2.57 

2.42 2.88 3.61 2.88 2.42 , (8) 

2.57 2.57 2.88 3.61 2.57 
3.03 2.57 2.42 2.57 3.61 y 

0.00 2.05 1.82 2.05 2.74^ 

2.05 0.00 2.51 2.05 2.05 

1.82 2.51 0.00 2.51 1.82 • (9) 

2.05 2.05 2.51 0.00 2.05 

2.74 2.05 1.82 2.05 0.00 



C/rjJ = 0.675) 



C/r„(J = 0.675) 



V' 
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Dynamical mean field (DMFT) calculation Avith frequency dependent interaction 



The Coulomb matrix elements ^ and (IS]), their frequency dependence (encoded by ImC/(w)) and the Hamiltonian 
matrix H^ for the p and d bands in the Wannier basis are the input of the DMFT calculation. DMFT neglects 
the momentum dependence of the self-energy and replaces the lattice problem by a 5-orbital impurity model with 
frequency dependent local interactions, and a self-consistency procedure involving Hk, which fixes the hybridization 
functions [S]. (How a momentum-dependent self-energy affects these multi-orbital calculations is an interesting open 
problem.) 

Since the LDA bandstructure [Hk] already captures some correlation effects in the d-orbitals, we modify the self- 
energy Yid by an orbital independent shift (double-counting correction) of 



Edc = i^°(nd - 1/2) - J{nd/2 - 1/2), 



(10) 



with rid the self-consistently computed number of d electrons. This double counting was found in Ref. 1121 to yield the 
best agreement between charge-selfconsistent and non-selfconsistent calculations. Note that the static value of the 
interaction appears in E^c ■ This is because the addition of bosonic modes in a dynamically screened model requires 
the introduction of additional double counting terms for these modes, which eliminate the bare interaction from the 
double counting formula. 

The new feature, compared to previous LDA+DMFT simulations, is the treatment of the full frequency dependence 
of the interaction. We employ the method developed in Refs. [T31 [TU which is based on the hybridization expansion 
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approach T5]. This technique can be easily generahzed to multi-orbital systems. We define the bosonic factor 
b{uj', t) = cosh [(r - f )w'] / sinh [^] and the function (vahd for < t < /3) 

Kir)^ r^^^^[5(.',r)-6(.',0)], (11) 

which (up to the sign) corresponds to the twice integrated nonlocal interaction. The frequency dependence of U then 
enters the hybridization expansion calculation in the form of a non-local interaction between each pair of creation 
and/or annihilation operators (irrespective of orbital) of the form Wij = cxp[siSjK{Ti — Tj)], where t,; > tj are the 
positions of the two hybridization events on the imaginary time interval and s = 1 for creation operators and —1 
for annihilation operators. Apart from these modifications, the impurity calculation proceeds as usual (with the 
static value of U for the evaluation of the interaction contribution to the Monte Carlo weight), via random insertions 
and removals of pairs of hybridization operators. For models with density-density interaction, this method is highly 
efficient. At the lowest temperature, T — 145 K, we used about 10 CPU hours per iteration. 

Analytical continuation 

The DMFT calculation yields an imaginary-time Green function which contains all the dynamic features encoded 
in the retarded interaction U{lu). However, to analyze the effects of the retarted interaction on the spectral function, 
an inversion procedure is required which allows us to go from the imaginary time to the real time domain. The 
stochastic noise in the CTQMC data makes this problem extremely difficult for standard Maximum entropy methods 
if one aims at resolving intermediate-to-high energy features like the satellites discussed in the main text. In Ref. fT6l 
a new analytical continuation procedure has been proposed which is based on the exact atomic limit properties of 
quantum impurity problems with retarded interactions. In that limit, the exact Green function becomes the product 
of a purely static (local in time) interacting part, Gstatic(''')j and the factor B{t) = Gxp[— _ftr(r)] with bosonic symmetry 
containing the full retarded tail of the interaction U. Gstatic and B are analytically known, and B is responsible for 
both the low-energy renormalization of the Green function, and the satellites resulting from the screening processeses 
embedded in U. The information encoded in the B factor gives the correct asymptotics and intermediate-to-high 
energy properties even away from the atomic limit, as the hybridization affects mainly the low-energy part of the 
spectral function. 

Motivated by these considerations we introduce an auxiliary Green function GauxiUary satisfying 

Gir) = Gauxiliary(r)B(r). (12) 

GauxiUary describcs mainly the low-energy features of the full Green function, thanks to the energy scale separation 
in the spectrum. Therefore, the standard maximum entropy method can be applied reliably to compute the spectral 
function Paux of GauxiUary i whilc pg - the spcctral representation of the Bose factor B - can be obtained via the 



numerical integration of Eq. ( 11 ) at any desired accuracy. The full spectral function is obtained from the integral 



1 _|_ g-/^" 
piiu) = / de pB{e) ^^ ^ g-^(,-^))(^ _ g-^e) Paux(a^ - e). (13) 
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The bosonic spectral function ps for BaFe2As2 is shown in Fig. lb of the main text. It essentially inherits the 



structures from IniC/(aj)/w . At zero temperature, Eq. (13) reduces to the convolution of ps and paux- If /Oaux('^) has 



a sharp peak at a; « 0, this convolution will produce a satellite at each of the energies corresponding to sharp features 
in pB- 

To analytically continue the self-energy, E(za;„), we define an effective Green function G^(icj„) = [iojn — Mcff — 
S(iaj„)]~^ and apply the above procedure to G^(ia;„). From its spectral function p^{uj), the calculation of S](w) 
follows straightforwardly through the Kramers-Kronig relations. 
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